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ABSTRACT 


In  the  present  paper  we  extend  in  two  wavs  some  results'  presents  in 
Kotz  and  Johnson  (Com.  i_n  Statistics  ( 1 982 ) ,  A1J_  }  relating  to  the  study 
of  distributional  aspects  of  effects  of  errors  in  inspection  sampling; 

(1)  Multistage  sampling  with  k  successive  samples  involving  the  possibil¬ 
ity  of  two  types  of  errors  in  inspection  (classifying  a  defective  individual 
as  non-defective,  or  a  non-defective  as  defective);  (2;  Single-staoe  sampling 
considering  several  types  of  defects  of  which  only  one  is  tested  on  inspection. 
Both  (1)  and  (2)  lead  to  novel  multivariate  distributions.  Their  structural 
properties  are  analysed  in  some  detail  and  some  applications,  in  particular 
those  in  Quality  control  are  discussed. 
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1.  INTRODUCTION 

We  have  recently  studied  effects  of  inaccuracies  in  inspection  on  the  proper¬ 
ties  of  acceptance  sarpling  procedures  (Johnson  &  Kotz  (1981),  Kotz  & 

Johnson  (1982a, b)).  In  particular  we  have  considered  two-stage  sampling, 
wherein  up  to  two  successive  samples  of  sizes  n^  ,n^  may  be  taken  fror  a  lot 
of  size  N,  containino  D  defective  items.  (The  second  sample  is  taken  only 
if  the  number  (Z,)  of  items  found  to  be  defective  is  within  certain  limits 
defined  by  the  sampling  schemes.)  It  was  supposed  that  inspection  is  not 
perfect,  so  that  some  defective  items  may  not  be  noticed  as  such,  whi 1  ■  some 
nondefective  items  may  be  classified  as  ’defective'. 

Similar  investigations  have  been  described  in  quality  control  literature  by 
Hoag  et  al.  (1975)  Dorris  and  Foote  (1978)  and  Rahali  and  Foote  (1982)  (see 
also  Armstrong  (1982)).  In  tiie  present  paper  we  shall  extend  this  work 
along  two  directions,  each  of  which  introduces  apparently  novel  multivariate 
distributions. 

First  we  shall  generalize  our  result  to  multistage  sarpling  with  k  successive 
samples  of  sizes  n^ ,  n^,  ...,  and  suppose  that  for  the  i-th  sample 

p^  =  probability  that  a  defective  item  is  correctly  classified 

p.j  =  probability  that  a  nondefective  item  is  classified  as  'defective'. 

Secondly,  returning  to  single-stage  sampling  we  shall  consider  the  case  of 
two  types  of  defects  one  of  which  is  relatively  easy  to  detect,  and  investi¬ 
gate  the  distribution  of  the  second  kind  of  defect  when  selection  is  on 
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the  basis  of  the  first  kind.  As  before  we  will  suppose  that  denotes  the 
probability  of  correct  classification  of  a  defective  item  and  p^  the  probabil¬ 
ity  of  classification  of  a  non-defective  item  as  'defective'. 

2.  MULTISTAGE  SAMPLING  ANALYSIS 


We  will  be  interested  in  the  joint  distribution  of  Z-j  ,Z^ , . . .  ,Zk ,  the  numbers 
of  items  classified  as  defective  as  a  result  of  inspection  of  the  sequence  or 
k  samples.  This  distribution  seems  to  be  new.  and  presents  some  features  of 
theoretical  interest,  as  well  as  having  the  possible  practical  applications 
we  have  indicated.  Chief  among  the  latter  is  the  possible  calculation  of 
acceptance  probabilities  for  multi-stage  (or  even  sequential)  sampling 
schemes  under  imperfect  inspection.  Some  practical  comments,  in  the  two-stage 
case  are  given  in  Kotz  &  Johnson  (1982b). 


Let  Y^,Y^,..:,Yk  denote  the  actual  numbers  of  defectives  in  the  1st,  2nd,...,k-th 
samples  respectively.  The  joint  distribution  of  Y  =  (Y^,...,Y^)  is  a  multi¬ 
variate  hypergeometri c  with 


m  i  k  n .  :N-::n.') 

(0-yi<ni(i=l  ,...,k);  D-N  +  ^  =  ]  n.  ■  Jk.  =  ]  y.). 
Symbolically,  we  write 


(1) 


Y  Mult.Hyppk(n;D,N) . 


We  note  that 

k  (r. ) 
E[  H  (Y.  1  ( 
i=l  1 

where  a^  = 


ni"Yi ) 


(s,) 


)]=  D 


(>:r4) 


(N-D) 


*i> 


k  (ri',si)  ( >-(  ri+S  -  ) ) 
I!  n.  1  1  }/N  11 

i  =  l  1 


a (a -1 ) . . . (a-b+1 )  is  the  b-th  descending  factorial  of  a. 


(2) 


Conditionally  on  Y,  the  Z's  are  mutually  independent,  with  Z..  distributed  as 


-3- 


the  sum  of  two  independent  binomial  variables  with  parameters  (Y.  ,p.)  and 
(n^-Y-.p!)  respectively. 

k  k  z.  y.  n.-y.  ,  z.-h  y.-h  n.-y.-z.+h 

pr[ini(zi=2i);W  =  i3i[vh:0(hM(z;v)pNi1  o-p,) 1  (l-o:)1  1  1  ] 

k 

=  r  b(z.;y.,pi;n  -y  p!)  (3) 

where  b(z^’,y^,p.;  n.-v., pi)  is  the  probability  function  for  the  convolution 
of  the  two  binomial  distributions  with  parameters  (y.,p^)  and  n.-y.,p^) 
respectively. 

The  unconditional  distribution  of  Z  is 


(N-  n.  n  k  ^n . ■  k 

-y  /'Id-:,':  »(Ji;yi.Pi;nryj,pp 


M  i  nn-  m  •  \  r. 

M2.J]  ■  (  v  <(  1  :  r 

J  l  -k1  ^ i  i  =  l  i  i-l 


The  limits  for  y  are  as  in  (1). 


Symbolically 


jBin(Y1,p1)*Bin(n1-Y1,pp 


|Bin(Yk,pk)*Bin(nk-Yk,p^) 


Mult.Hypg  (n:D,.\') 


where  *  stands  for  convolution  and  is  the  compounding  operator  (as  defined, 
for  example,  in  Johnson  &  Kotz  (1969,  p.  18A)). 

This  distribution  might  be  called  "Multivariate  Hyperqeometric-Convcluted 

Binomial (s) 


The  conditional  (on  Y)  r^-th  factorial  moment  of  Zi  is 


(r«)  r.  r.  .  r.-h  ,  (r.-h) 

£tz(  ’  in  •  5h:0(h’)pbp;  ’  v(h>(Vv  1  . 
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The  unconditional  r  =  (r . , . . . ,r^) -th  factorial  moment  of  Z  is 

k  (r.)  k  (r. )  k  (r.) 

Ur(2)=  E[  n  z.  1  ]  =  E  E[  r  z  1  |Y]  =  £  [  ::  E[z  iy]] 

X  i=l  1  i  i=l  1  -  i=l  1 


>.-h) 


k  r.  r.  .  r.  -h 

■  Ek  '  ’P  Pi  ’  ‘VV  ° 
1  =  1 


(6: 


=  V 


k  (r. ) 

1  r.  rb  k  r.  ■  .  r.-.  .  ('  .  )  (  (r.-\)) 

1=1  1  k  /  1  W ,  1.  .1  1  r,  1  /„  1  1 


(:r.) 


=0 


i '  1  -  k  i  =  l  i 


(  ')[•:  V’  ’  &  1  (*:-D) 

_•  i  ■  1  1 


k  (r. ) 

\-,ni  r  ()  (V) 

•  ^rtr- 1.4  o< •»(»-» 

N  1 


r . 


16') 


where  g.  is  the  coefficient  of  x’  in  (ri-tn.x)  1  .  (Note  that 

i  *  1  1  1 

(i)  if  sampling  were  with  replacement,  (b)  would  be 


k  (rJ  k 
>lr](7)  =  '  :  n  1  [  dp  4  (n-p)p;  3/N 

w  i=l  1  i=l  1  1 


r.  .  r. 

.  h/»,  1 


(ii)  if  p.  =  p  and  p.!  =  p'  for  all  i -correspondinq  to  constant  quality  of 
inspection  throuohout  -  then 

9.  =(  ' i)p'p‘  ^ 


In  particular 

E[Zi  ]  =  ni[DN'1P14(l-DN'1)pi'  1  -  n.'P]  (7.1) 

where  =  DN_^pl- +(1  -  ON -  ^ )  p  1  is  the  probability  that  an  individual  chosen  at 
random  in  the  i - th  sample  will  be  classified  as  defective  (whether  it  really 
is  so,  or  not).  Also 
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var(Zi)  =  ^p.d-p.)  -  n.  (n.-l )  (N-1  rV  ( 1  -DfT1 )  (p.-p!  )2  (7.2) 

and 

cov(Z^ ,Zj )  =  -n1n.(N-1)'1DN'1(l-DM‘1)(pi-p!)(pj-p:)  (8) 


If  p.  >  p.'  and  p.  >  p‘.  as  one  would  hope,  if  inspection  is  to  be  of  any  use 
^  *  JO 

at  all,  the  covariance  is  negative,  as  might  be  expected,  since  the  covariance 

of  Y.  and  Y.  ’s  negative.  If  p.  =  p!  or  p.  =  u'  so  that  it  is  irrelevant 

i  J  *  0  0 

whether  an  i ten  in  the  corresponding  sa:'iple(s)  is  defective  or  not,  then  the 
covariance  is  zero  -  in  fact  the  correcoondinn  Z  (or  Z's)  is  independent  of  all 
other  Z's. 

It  is  easy  to  write  down  the  conditional  distribution  of  Z.,  given  Z.,  but  it 

i  J 

is  rather  complicated  in  form.  We  can,  however,  derive  the  regression  function 
in  the  fol lowing  way. 


We  have 

EtV'V  -  Vi  ♦  <vVi 

so 

EtZ^Zj]  ■  P^iZj]  *  I-itn^Etg  ZjJ) 

(9) 

Now 

E[Vi  :YjJ  ■  n^D-rl/fN-np 

so 

EOpZj]  ■  n.lN-^rho-EfYjiZ^) 

(10) 

In  order  to  evaluate  E[Y..  ;Z.]  note  that,  for  the  j-th  sample, 

u  J 


and 


Pr[item  is  defective j item  classified  as  defective]  =  DN  p./p. 

J  J 


Pr[item  is  defective! item  classified  as  nondefective]  =  DN-1 ( 1 -p . )/( 1 -p . ) 

J  0 


SO  ECYjlZj]  =  DN_1[Zj(Pj/Pj)  ♦  (nj-ZjJIO-PjJ/d-Pj)  *3 


(ID 


From  (9),  (10)  and  (11)  we  obtain 
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iUiUj]  --  niP!  +  (prp!)  ECY.jZj] 

=  n.p^  +  ni(pi-p!)(N-nj)‘1(D-E[Yi  Z ^ ] ) 

i  (1-Drr1)(n.-P!) 

=  n.p!  +  n  -  ( P- -p-  )DN  1  1 - J — J- 

(N-n^d-fK) 


The  regression  is  linear,  and  the  siqn  of  the  regression  coefficient  is 
opposite  to  that  of  (p.-p ! )  (Pj'Pj )  anreeir.q  with  the  sign  of  cov(Z^Z^)  in  '  P, ) 


3.  MULTISTAGE  SAMPLING:  A  SPECIAL  CASE 


If  we  take  n^  =  n^=  ...  =nk  =  1  we  have  the  first  k  stages  of  a  fully  sequent¬ 
ial  sampling  procedure.  In  this  situation  the  only  possible  values  of  each  of 
the  Y‘s  and  Z 1 s  are  0  and  1.  Formula  (3)  becomes 

NT  M  k  k  ]-Zi  zi  ]'zi 

P^3  "  Q  vy  •••/»(?-  v  >  .''Vi  <V>  Mi-y^p:  (1-p;)  03 

J 1  -  k  -  v- 1 


Collecting  together  terms  with  the  same  value  of  y.  =  y  (corresponding  to  the 
total  number  of  defective  items  selected)  we  net 

mm]  ■  (J)'1  :kC>  ir,  ,k  p/pipi-Pii'd-p;) f  o; 

~  -  i=l 

where  summation  with  respect  to  _i  =  (  . k)  and  ^  ^  ^ )  is 

constrained  by 

.(  /.+;•.)  =  y,  and  +o!)  =  ,z.  =  z 

'i  =  l'  l  i /  i  -- 1 '  i  i  '  i  =  1  l 

(corresponding  to  the  total  number  of  items  classified  as  defectives)  and  for 

each  i,  one  of  i  ,i; ! )  is  1,  the  other  three  are  each  zero. 


Yet  another  way  of  writing  (13)  is 


-7- 


Pr[Z=z]  =  (p)'1  ry(o"y) ‘[coefficient  of  xZuV  in 
k 

n  ;p.XU+P-X+(l-p,)u+l'P^  ]  (13)" 

j=l  J  J  J  J 

It  may  be  noted  that  this  shows  that  Z.  is  a  sufficient  statistic 
(for  D,  supposing  N,  p's  and  p"s  are  known).  ]f  p.  =  p  and  p!  -  p'  for 

7  v 

all  i,  we  have  the  coefficient  of  x  u-  in  pxi.  +  p'x+( l-p)u+(l-p‘ )  •  or  tne 
right  hand  side  of  (13)'*  . 

4.  ASSOCIATED  DEFECTS:  ANALYSIS 

We  now  suppose  that  there  are  two  types  of  defects  -  (I)  and  (2)  -  ana  that  in 
a  population  of  size  N  there  are  Dq^  individuals  with  g  type  (1),  and  h  type 
(2)  defects  (q,h  =  0,1).  (Of  course  D^+D^+D^+D^  1  =  N.)  For  example  (1) 
might  represent  surface  irregularity,  with  (2)  cerrespondi nq  to  internal  flaws 
in  material,  such  as  metal  bars  or  plates.  In  many  situations  the  different 
types  of  defect  correspond  to  different  modes  of  failure. 

A  random  sample  of  size  n  is  taken  (without  replacement)  from  the  population; 
and  each  of  the  chosen  individuals  is  examined  for  presence  of  type  (1) 
defect.  We  are  interested  in  the  number  (Z^)  of  individuals  classi f ied  as 
possessing  defect  (1);  and  among  these  Z^  individuals,  the  number  (Z^) 
actual ly  possessing  defect  (1)  and  the  number  (7p  possessing  defect  (2). 

The  distribution  of  Z^  is  given  in  Section  2,  and  also  in  Kotz  &  Johnson  (1982b). 
The  distribution  of  Z*  is  the  hypergeometric-binomial 

Binomial  (Y.p^ )  A  Hypq.  (,4) 

The  probability  that  the  sample  will  contain  yQh  individuals  with  g  type  (1) 
defects  and  h  type  (?)  defects  is 
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1  1 


1  1 


D 


Pr£  0  ■  (vQh=v  h)]  =  [  ’  r.  (  Qhn/Q  (y00+y0i+yio+yn  =  n) 

g=0  h=0  9  9  g-0  h=0  ygh  n  Ob  Ul  IU  II 

(A  multivariate  hypergeometri c  (n;D,N)  distribution  -cf  (4a)  where  n  is  a 
vector  but  D  is  not.)  For  the  -  ^l00"l01’’10  joint  factorial  moment 

of  Y  we  have  the  eApresssion 


(15) 


/  \  i  i  !  QJ  (  \ 

u,  A v)  =  n  '>[  :  D  .  ]/hl  '  where 

'  q=0  h=0  9 


'00+'Cl+  io+-'n  {16) 


Given  y  r  ^y00’'V01  'v10*yl  1  ^  ’  Z1  1S  distributed  as  the  sum  of  two  independent 
binomial  variables,  with  parameters  vy-|Q+y-] -|  .P-] )  and  ^y00+y01,pi^  r®specti  vely . 
The  first  of  these  variables  is,  in  fact  Z^.  Similarly,  Z*  is  distributed  as 
the  sum  of  two  independent  binomial  variables  with  parameters  (y^,p.)  and 
(y01’pl^  resPecti vely.  Introducing  four  independent  binomial  variables  W  h 
with  parameters  (y^.  :  +  (l-:^)pj)  for  q,  h=0,l  with  vq  =  i  1  - ( - T ) Q  .  ,  then 

conditionally  on  Y  =  ^ 


71  =  W00+W01+W10+Wll 

(17.11 

Zt  =  W10+Wll 

(17.2) 

Z*  =  WQ1+Wn 

(17.3) 

We  have 

E[Zliy]  =  (yi0+yll)Pl+(y00+y01)pl;var{Zliy)=(y10+yll^pl(1 

-pi)+(yoi+yoo)pi(1- 

(18.1) 

Etz^ly]  =  y1iP1+y0ipi;var(z2iy)=ynpi(1'pi)+yoipi(1'pi#) 

(18-' 

cov(Z^  »Z||y)  =  var(W1Q+W1 1 |y)  =  var(Z£]y) 

(18.3) 

corr(Z^  »Z^|y)  =  (var(Z||y)/var(Z1  [y)  /•' 
cov(Z^,Z|jy)  =  var(Wn)y)  =  ynP](l-p1) 

(18.4) 

(18.5) 

(18.5) 
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corr(Z^,Z|!y)  = 


r  2  v 

I  ynpi^"pi^  ^  _  i 

i  ^1  o+y  n )  -y  l  i  pi  ^  “pi  ^+yoi  PT^  "pv^ 1 ! 


(18.6) 


The  distribution  of  Z£  is  the  multivariate  hyperqeometric-convol uted  binomial 


!>)  >  /\ 
i}  Yir 

Moments  can  be  obtained  similarly  as  in  Section  2. 


(Binomial  (Y^  ) *3i nomi al  (y10>Pj )  y  XjY  Mult.Hypg(n;D^  ,01Q;N)  (19) 


For  an  individual  chosen  at  random  from  the  population 

Pr[classified  as  having  defect  (1)]  =  N*1  •'  (Dgi  +  DQ0)pj  +  (D10+D1 1  )p^ 

=  n'1  (D0ni+Di  P])  r  P] »  sa-v  (i9.i : 

(D^  =D.jQ+D.ji  =  total  number  of  individuals  in  the  population  with  defect  (1); 

Dq  -  N-D,  ) 


and  , 

^  ( 9qi o. +Di i Pi ) 

Pr[havinq  defect  (2)  classified  as  having  defect  (1)]  =  - - - - 

pl 

P2  1 

=  —  ,  say  (19.2) 
P1 

Hence  ECZ^Z-j]  =  Z]  P2  i /Pi  (20.1) 

and  (from  (7.1)) 

E[Z*}  =  np2,{  (20.2) 

Formula  (20.1)  is  also  valid  conditionally  on  Y=y,  so  we  have 

covUj.Z*)  =  ^  E[zf ]  -  n26P2jl  (21.1) 

and  from  (7.2)  and  (7.1 ) 


np?  |  -I  n  1  ^1  ? 

cov(ZrZ*)  =  -Hi  -^(1-^)  -  g-f  4-(1-  V.)fPrPj)  ■  (21-2 

pl 

We  must  have  cov(Z-|,Z|)  >  0  because  from  (20.1),  F[Z£  Zj]  is  an  increasing 
function  of  Z-j . 


Extension  to  situations  in  which  there  are  m  (  1)  types  of  defects  -  (2), 

(3),  (m+1)  -  in  addition  to  the  type  (1)  which  is  inspected  direct'.,  s 

stra  ightforv/ard.  The  only  essentially  new  problems  are  the  joint  distribute 
of  -  the  numbers  of  individuals  with  type  (2) ...  .,(;•  +  ! )  defects 

among  those  classified  as  having  defect  (1);  and  also  the  distributions  o'- 

variables  like  the  number  anong  these  individuals,  with  both  (i)  &  (j; 

type  defects  (i  /  j  .-2).  Using  an  obvious  notation  (with  subscripts  0(1) 
indicating  absence  (presence)  of  the  corresponding  type  of  defect)  we  have, 
for  example,  corresponding  to  (19.2) 


Pr[having  defects  (2)  and  (3),  but  not  (4),. ..(m+1)  classified  as  having 


defect  (1 ) ]  = 


Dono...opi  *  pi iio. .^o_pi 

°0. . .pi  +  Dl...  P1 


(22) 


where  D_ 


a^O 


vi*°  t**-* 


m+1 


(cj  =  0,1)  (D  is  the  quantity 


previously  represented  by  P  . ) 

Considering,  for  simplicity,  the  case  m  =  2,  we  now  obtain  an  expression  for 
cov(Z|,Z^).  We  have,  analonously  tu  (17.3) 


z2  "  wi 1 1  +  won+  wno  +  woio 
Z3  =  wni  +  won+  wioi  +  woci 


(23) 


where  the  W's  are  independent  binomial  variables  and  the  parameters  of  W 
are  frghf  Vl  *  n"’g)pi)- 


ghi 


Symbolically 
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IBin(yno,pi )  *B  1  n  (  yqi  o  *pi ) ' 


(24) 


2*\  i 

21 Bin.  (V1  ^  ,d1  )*Bin(YQ11  ,nj)-  \  i’/.Mult. Hype (n 


i  B i n ( V i  oi  >P-j  ;  ""B l  n  (  ^q01  ,P1  ^  ■  " 


(Here  Y  =  (Vin,  YC11,  Y110,  Yq10,  t]Q1.  y001  )  •  ! 

(s ?)  (s3) 

A  general  expression  for  the  joint  factorial  norent  E [ 1\  ]  car  be 


obtained  in  the  following  way. 

T(4) 


S  o )  ( s  o )  i  ,r  \  S  y  ( i  i )  ( i,  )  ( i  -3 )  ( i  /i ) 


Z2  ‘  Z3 


-  H 


s-,  (j,)  (j,)  (j,)  (jj 


2,,/J3/  _  ,r(C)/'2x11',l/11v,2V,3/I,v,4/.  -(4)  /11lJ2/l..'J3/1  ,w4 


on  "no  "oio  :  ,  (j  )wm  Sii  wioT  Soi 


„tA\  IA\  s?  ^1  J,  (jr<j)  (inn)  J2  (j2-h)  (i2+h) 

=  ‘w  out  (  1  n,  1  u  1  .  (  ^  ‘  u  ‘ 

n  i\t  -  _  in  wm  1h  ' 1 2  won 

n=0 


1  J 


'I  "i  a=0  0  1 


do)  (i4)  (jo)  (jJ 

+  u  J  ll  J  W  ^ 

"no  oio  "ioi  ooi 


(25; 


where 


•(4)  .  - 


i ,  i  „  i ,  i 
1  c  a  ■ 


_  (4  \  . 

J'  -  _  and  utilizing  the 

j  J  i  3  2  -1 3  j  4 


Vi 


U)  U)  (,,)  U)  _‘2  (1?-U)  (,,+u) 

identity  a  a  =  a  (a-.»,+  ,)  --  (  ;  ,  a 

1  1  u=0  u  1 


(s2)  (s,); 


An  expression  for  E[Z:$  Z^  ;Y-y]  is  obtained  by  replacing  W^(  in  (25) 

by  Yg^]Pg  where  p^  =  g*P]  +  (l-'q)nj.  Then  taking  expectations  with 
respect  to  Y,  we  obtain 

1  2  s„  So  j,  jo  9+ii+io+Jo  h+i o+jo+j «  (26) 

1  V  I  £\t  -1  i  /  •  w  I  d  d  .  d  d  4 


(So)  (So)  , (4\f4x*l  “2  s?  S3  j,  jo  g 

e[z*  2  z*  3  ]  -  y(4,n4)i  y  (i2)(i3H.1)(  2)Pl 

£  i  i  j  g=o  h=0  ~  i  q  n 


'1 


(jrg)  (j2-h)  n(s2+s3+g+h)  ( i -j  +g )  (i2+h)  (i3)  (i4)  (j3)  (j4) 

(so+So+o+hJ  Dlll  D01 1  D1 10  D010  D101  D001 


*  d  ’2 


If  m  >  2,  Dgh-  is  replaced  by  D^hi  ;  joint  factorial  moments  of  three  or 
more  Zf's  (i  >  1)  can  be  obtained  by  similar  techniques,  though  the  formulas 
rapidly  become  more  cumbersome. 
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For  lower  order  moments,  direct  calculation  is  often  simpler  than  using  general 
formulas . 


We  now  outline  the  calculation  of  cov(Z^.Z^).  We  have  from  (23)  and  (24) 


cov(7-.z*^)=  y111P1n-P1)+y011Pid-p] ) 


( 27 ) 


whence 


Etz2z3  ympi^"pi  '+yoiipi(1'pi )+  (ym+yi io^pi  +  "yoi i^yoio^Di 


so  (using  (15)) 
cov(Z|,Z*)=E[2JZ|]-E[2|]E[2*] 


vym+yioi  )ri+(yoir-vooi)pi  ’ 


n^di n pi ( i ~pi ^+Doi i pi ( 1  _pi )i+  In  '  n(nT)  (Dmpi+Donpi  ' 


+{S  IF-TT  "  J"2  ■'(Dni+Dno)pi  +  (Don"DoioDi)  :'  (Dm+Dioi)pi+(Don+Dooi)pi 


n-l 


n{diiipi(1‘  n-i  pi)+Donpi(1'  n- 


n-l  , 

M  ■]  P1 


(23) 


"  n2(n-i) i(Di11+Dllo)Pl  +  !D°ii+Do1o)|,'i  (Dm+Dioi)pi+(Don+Dooi  M 

This  covariance  can  be  positive  or  negative.  Bo.h  l*  and  Z^  are  positively 
correlated  with  Z-j  (the  total  number  of  individuals,  classified  as  ’defective', 
of  which  Z|  and  Z^  are  subsets),  but  they  may  be  negatively  correlated  in  the 
population.  The  latter  situation  corresponds  to  values  of  0qlQ  and  D  ^  which 
are  large  relative  to  D  ^ .  (g  =  0,1). 

When  N  is  large  compared  with  n, 

cov(Z£,Z*)  =  n[p1 ; t t ~ ( P1  ;n+  Pj  lo^hlll*  pl|01^' 
where  p,  (M  *  N',(D,htP,*°ohtpiK 
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Taking  as  an  example  P1  =  0.90,  pj  =  0.10 

Dni/N  =  0.01;  Dno/N  =  0.1  ;  D^/N  =  0.15 
D0,,/H  =  0.01;  D0,0'N  =  0.15;  0Q01 /N  -  0.0 

we  have 

p-j  | <1 1  =  0.010;  p^  i  =  0.105;  ^  Q]  =  0.155 

so  h  n"(pi  lTh  10^1  n+pi  oi^  0 

corresponding  to  negative  cov(Z^,Z*)  . 

On  the  other  hand,  if  D^^/N  and  Dq^/N  are  each  increased  to  0.05,  the  other 
parameters  remaining  the  same,  we  have  ^  =  0.050,  while  p^  ^  and  p^ 
remain  unchanged,  so 

Pi  ]ii_  ( P]  ii+  P-|  ;-|q)  (P-|  i-|i  +  P |  Q-j )  -  0.050  -  (0.^55  x  0.205)  .  0 
corresponding  to  positive  cov(Z£,Z*). 

The  same  formulas  apply  when  there  are  m(  2)  types  of  defect  other  than  (1), 
replacing  Dghi  by 

At  the  cost  of  some  elaboration  in  the  formulas,  we  can  allow  for  the  possibi¬ 
lity  that  presence  or  absence  of  a  defect  of  type  (2)  may  affect  the  probabil¬ 
ity  of  correct  classification  in  regard  to  defects  of  type  (1).  Introducing 
the  notation 

2^1  ^1  ^  f°r  Pr°bat>1 1  iti es  of  detection  of  (1)  in  the  presence  (absence) 
of  (2) 

and 

2P-j  for  probabilities  of  incorrect  assignment  of  (1)  when  no  (1) 

is  present.  In  the  presence  (absence)  of  (2)  we  would  still  have  a  model 
of  form  (17)  but  the  parameters  of  the  binomial  distributions  of  the 

VS  W°Uld  n°W  be  *yg0’V2Ml  +  (1"V2pi'  for  h  -  0;  g  =  1 ,2 
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g*  2P1  +  (1"q)5pi)  for  h  =  1 ;  g  =  1 ,2  . 
The  probability  of  having  defect  (2),  if  classified  as  having  (1)  would  be 


00l’2pi  +  Dll‘2Pl 


p2  1 


3oo'?pi  +  loi*2pi  +  3io'?pi  +  Dn*2pi  p* 


(30) 


where  p*  =  ^(Dgg-^pj  +D01*2pi  +  D10'2P1  +  D11  *2c'l  ^  =  Pr[classified  as  havin' 
defect  (1 )] 

and  p*;1  *  ^Doi'2pl  +  Dll‘2plJ 

Formulas  (20.1)  and  (20.2)  would  still  be  valid,  with  ^  replaced  by 

p*,^  respectively. 


Pi. 


Extension  to  situations  with  m( --1 )  types  of  defect,  other  than  the  one  (type 
(1))  which  is  inspected  directly  is,  again,  straightforward. 


In  view  of  the  model  (17)  which  applies,  with  appropriate  adjustments ,  to  all 
the  cases  mentioned  above,  the  joint  distribution  of  Z's  and  Z*'s  is 
asymptotically  multinomial  as  the  population  size  N  increases  indefinitely, 
with  the  ratios  n:D's:N  remaining  constant,  or  tending  to  fixed  values. 


5.  SOME  APPLICATIONS 


Although  this  paper  is  concerned  primarily  with  some  novel  compound  multivariate 
discrete  distribution  which  can  arise  in  connection  witn  faulty  inspection 
rather  than  in  specific  applications  ,  we  shall  indicate  in  this  section  a  few 
circumstances  in  which  knowledge  of  these  distributions  may  be  useful  and 
directly  applicable  to  specific  investigation'  and  inquiries. 

The  results  in  Section  2  are  relevant  to  studies  cf  robustness  of  multistage 
Simple  procedures  to  errors  incurred  in  inspection  and  consequently  to  the 


actual  construction  of  such  procedures.  They  would  also  be  directly  relevant 
to  construction  of  tests  for  inequalities  anonq  the  p.'s  and/or  p!  's  which 
could  be  one  aspect  of  attempting  to  detect  the  existence  of  faulty  items. 
(Evidently  if  p^  and/or  vary  with  i,  they  cannot  be  identically  equal  to  1 
or  0  respectively  for  all  i).  The  distributions  derived  in  this  paper  are 
also  indirectly  relevant  to  construction  of  tests  of  hypothesis  of  no  faults 
(p  =  1,  p'  =0)  assuming  p.,  p.!  do  not  depend  on  i.  Some  attempts  in  this 
direction  have  been  made  in  Johnson  &  Kotz  (1982)  while  analyses  of  ways  in 
which  cost  consideration  can  be  allowed  for  in  faulty  inspection  problems  are 
sketched  in  Kotz  &  Johnson  (1982b). 

The  results  in  Section  4  are  relevant  to  assessment  of  performance  of  proce¬ 
dures  for  identifying  individuals  with  defects  of  type  (2), say  (especially 
in  those  cases  when  these  defects  are  not  easily  detectable)  by 
observing  the  existence  or  non-existence  of  defects  of  type  (1),  and  the 
robustness  of  this  assessment  to  actual  numbers  of  faults  among  inspected  iters. 
In  these  circumstances  it  may  sometimes  be  appropriate  to  carry  out  a  IOC. 
inspection  -  that  is,  n  =  N  -  though  the  more  general  formulas  we  have 
derived  are  of  course  of  greater  flexibility  and  are  useful  in  various  situa¬ 
tions  when  total  inspection  is  either  not  feasible  or  too  costly.  Indeed, 
studies  in  this  direction  will  involve  introduction  of  cost  functions  allow¬ 
ing  for  costs  of  sampling  and  losses  due  to  erroneous  retention  of  defective 
individuals  of  type  (2)  and  the  erroneous  rejection  of  non-defective  type  (2) 
items.  See  Kotz  &  Johnson  (1982b)  for  an  appropriate  model  and  some 
preliminary  results  for  the  case  of  two-stage  sampling  with  defects  of  a 
single  type.  Finally  questions  of  choice  of  which  type(s)  of  defects  to 
inspect  for  can  also  arise  in  this  context. 
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